Rapid sympathetic cooling to Fermi degeneracy on a chip 
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Neutral fermions present new opportunities for 
testing many-body condensed matter systems, re- 
alizing precision atom interferometry, producing 
ultra-cold molecules, and investigating fundamen- 
tal forces. However, since their first observation 
Pi, quantum degenerate Fermi gases (DFGs) have 
continued to be challenging to produce, and have 
been realized in o nly a handful of laboratories 

USSHIESSSEi- In this Letter, we report 
the production of a DFG using a simple appa- 
ratus based on a microfabricated magnetic trap. 
Similar approaches applied to B ose -Einstein Con- 
densation (BEC) of 87 Rb [T3, E3 have acceler- 
ated evaporative cooling and eliminated the need 
for multiple vacuum chambers. We demonstrate 
sympathetic cooling for the first time in a micro- 
trap, and cool 40 K to Fermi degeneracy in just 
six seconds — faster than has been possible in 
conventional magnetic traps. To understand our 
sympathetic cooling trajectory, we measure the 
temperature dependence of the 40 K- 87 Rb cross- 
section and observe its Ramsauer-Townsend re- 
duction. 

Microfabricating the electromagnets used to trap ultra- 
cold atoms leads to a series of experimental benefits. 
Decreasing the radius R of a surface-mounted wire in- 
creases the maximum magnetic field gradient as Br 1 ! 2 
|l3j |. Since the oscillation frequency u> of the trapped 
atoms increases linearly with transverse field gradient, 
decreasing R from centimeters to micrometers can in- 
crease the confinement frequency by orders of magni- 
tude. In addition, one can envision a "lab on a chip," 
in which multiple devices are integrated on a single de- 
vice, expediting applications for complex manipulation 
of fermionic atoms for simulations of strongly corre- 
lated systems, quantum transport experiments, collision- 
insensitive clocks, and precise interferometry 0, Il5| . 
The strong confinement provided by a microfabricated 
electromagnet (/iEM) trap also has a practical advantage: 
it facilitates faster cooling, which relaxes constraints on 
vacuum quality and leads to a tremendous simplification 
over traditional DFG experiments that require multi- 
ple ovens, Zeeman slowers, or two magneto-optical traps 
(MOTs). 

In our system |16|. the entire experimental cycle 
takes place in a single vapour cell (Fig. QJi). Counter- 
propagating laser beams collect, cool, and trap 2 x 10 7 
40 K and 10 9 87 Rb atoms in a MOT. Atoms are trans- 
ferred to a purely magnetic trap formed by external 
quadrupole coils and transported to the chip 5 cm away. 
Figure \I]p shows several microscopic gold wires sup- 



ported by the substrate. In the presence of uniform 
magnetic fields, current flowing through the central 'Z'- 
shaped wire creates a magnetic field minimum above the 
chip. At the centre of this trap, the 40 K radial (lon- 
gitudinal) oscillation frequency is u>±/2tt — 826 ± 7 Hz 
(w £ /2tt = 46.2 ± 0.7Hz). The corresponding 87 Rb trap 
frequencies are a factor of \J mub/rriK ~ 1-47 smaller, 
where mjy, and ttik are the atomic masses of 87 Rb and 
40 K, respectively. 

After loading, the 1.1-mK-deep chip trap holds ap- 
proximately 2 x 10 5 40 K and 2 x 10 7 87 Rb doubly spin- 
polarized atoms, at a temperature > 300 /iK. Lower tem- 
peratures are achieved by forced evaporative cooling of 
87 Rb. A transverse magnetic field oscillating at RF fre- 
quency ^rf (typically swept from 30 MHz to 3.61 MHz) 
selectively removes the highest energy 87 Rb atoms by 
driving spin-flip transitions to untrapped states. The 40 K 
atoms, with smaller Zeeman splittings, are not ejected 
but are sympathetically cooled 0, 0, 0] by thermal- 
izing with the 87 Rb reservoir via 40 K- 87 Rb collisions 

slam. 

The evolution of temperature T and atom number 
N during sympathetic cooling is measured by releasing 
atoms from the trap and observing their expansion with 
absorptive imaging. Figure [21 shows the cooling of 40 K 
and 87 Rb to quantum degeneracy. In the degenerate 
regime, bosons accumulate in the ground state (form- 
ing a BEC), while fermions fill the lowest energy levels 
of the trap with near-unity occupation. Fermi degener- 
acy can be quantified with the fugacity Z: the ground 
state has occupation Z / (1 + Z), which approaches 1 in 
the high Z degenerate limit and Z in the non-degenerate 
limit. Due to the tight confinement of the /iEM trap, 
cooling increases the 40 K fugacity by 10 12 in only 6 s. 
The steep ascent of 40 K in Fig. [3 also demonstrates the 
efficiency of sympathetic cooling. Since 40 K is a rare iso- 
tope and therefore more difficult to collect from vapour 
than 87 Rb, this inherent efficiency is a significant advan- 
tage over direct evaporation. To our knowledge, this is 
the first observation of sympathetic cooling, of Fermi de- 
generacy, and of dual degeneracy in a ^EM trap. 

Below T « 1 /iK, we observe two independent sig- 
natures of Fermi degeneracy. First, we compare the 
rms cloud size of 40 K and 87 Rb (or its non-condensed 
fraction) by fitting the density profiles to a Gaussian 
profile. As described in Methods, this is an appropri- 
ate method for finding the temperature of a classical 
Boltzmann gas. Figure [21 shows that the apparent (i.e., 
Gaussian-estimated) 40 K temperature approaches a fi- 
nite value while the 87 Rb temperature approaches zero, 
even though the two gases are in good thermal con- 
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FIG. 1: A simple apparatus for Fermi degeneracy, (a) The dual-species MOT (red sphere) is formed at the intersection 
of six laser beams. The cloud is then magnetically trapped using external quadrupole coils (blue) , transported 5 cm vertically 
using an offset coil (purple), and compressed in the fiEM trap, (b) Schematic diagram of the central region of the fiEM 
chip. A magnetic trap is formed 180 /im above the surface at the location marked with a white "X" by applying 1% — 2.0 A, 
Itj = 30 mA, Bbias = 21.4 G, and -Bioffe = 5.2 G. Wire widths from left to right are 20 /im, 60 /im, and 420 /im. 



tact. In fact, this deviation is evidence of the "Pauli 
pressure" expected of a gas obeying Fermi statistics 0: 
at zero temperature, fermions fill all available states up 
to the Fermi energy E-p = h^Nto^LOi) 1 ^ , where N is 
the number of fermions. For our typical parameters, 
E-p ps fcex 1.1 fxK. We plot data with thermal (dia- 
monds) and Bose-condensed (circles) 87 Rb separately, to 
show that the density-dependent attractive interaction 
between 40 K and 87 Rb does not affect the release energy 
significantly. A second signature of Fermi statistics is ev- 
ident in the shape of the cloud. Figure compares the 
residuals of a Gaussian fit (which assumes Boltzmann 
statistics) to the residuals of a fit which assumes Fermi- 
Dirac statistics. The Fermi distribution describes the 
data well, with a x 2 three times lower than the Gaus- 
sian fit. After all 87 Rb atoms have been evaporated, we 
use Fermi-Dirac fits to measure temperature, and find 
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FIG. 2: Sympathetic cooling in a chip trap. Spin- 
polarized fermions without a bosonic bath cannot be success- 
fully evaporated (blue diamond). However, if bosonic 87 Rb 
(red squares) is evaporatively cooled, then fermionic 40 K is 
sympathetically cooled (blue dots) to quantum degeneracy 
(greyed area). For bosonic 87 Rb, the vertical axis is the occu- 
pation of the ground state; for fermionic 40 K, the vertical axis 
is the fugacity, as discussed in the text. These two quantities 
are equivalent in the non-degenerate limit. A typical run-to- 
run spread in atom number is shown on the right-most point; 
all vertical error bars are smaller than the marker size. 



k^T/Ep as low as 0.09 ± 0.05 with as many as 4 x 10 4 
40 K atoms. 

We empirically optimize the sympathetic cooling tra- 
jectory, and find that RF sweep times faster than 6 s are 
not successful, whereas 87 Rb alone can be cooled to de- 
generacy m 2 s. This indicates that 40 K and 87 Rb rether- 
malise more slowly than 87 Rb with itself. Measuring the 
temperature ratio during sympathetic cooling (Fig. 
reveals that 40 K lags behind 87 Rb at high temperatures, 
despite the fact that our optimal frequency ramp starts 
slowly (when the atoms are hottest), and accelerates at 
lower temperatures. 

In the low-temperature limit, we do not expect the 
cross-species thermalisation to lag the 87 Rb- 87 Rb ther- 
malisation, since the 40 K- 87 Rb cross-section ORRb = 
1480 ± 70 ran 2 [H exceeds the 87 Rb- 87 Rb cross-section, 
0RbRb = 689.6 ± 0.3 nm 2 |2l|. However, several con- 
flicting values for <7KRb have been recently presented 

us Etna mill. 

To investigate CKRb further, we measure the cross- 
species thermalisation rate [24| at several temperatures. 
Starting from equilibrium, we abruptly cool 87 Rb by re- 
ducing i'rf, wait for a variable hold time to allow cross- 
thermalisation, and then measure the 40 K temperature, 
as shown in Fig.^p. We repeat this measurement at sev- 
eral temperatures, and fit each to the model of Ref. [25| . 
We find that the cross-section has a dramatic dependence 
on temperature (see Fig. 010, decreasing over an order of 
magnitude between 10 and 200 /iK. 

The simplest model for atom-atom scattering uses a 
delta-function contact potential. Fig. ^J: shows that the 
s-wave scattering cross-section of this "naive" model (fur- 
ther described in Methods) would predict a higher <7KRb 
than ORbRb throughout our cooling cycle, in stark con- 
trast to our data. Better agreement is given by an effec- 
tive range model |26j| , which includes a reduction in scat- 
tering phase (and thus cross-section) below the naive ex- 
pectation. Our highest temperature data point lies below 
the effective range prediction, however a more sophisti- 
cated analysis may be required to extract a quantitative 
measurement for this point, due to severe trap anhar- 
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FIG. 3: Observation of Fermi statistics. Due to Pauli 

pressure, Fermi degenerate 40 K clouds appear to stop get- 
ting colder, even when the reservoir temperature approaches 
zero. The apparent temperature of the fermions, as measured 
by Gaussian fits to images of 40 K clouds, is plotted versus 
temperature of both thermal (diamonds) and Bose-condensed 
(circles) 87 Rb. Data is compared to a Gaussian fit of theo- 
retically generated ideal Fermi distribution (solid line) and 
its classical expectation (dashed line). Both temperatures 
are scaled by the Fermi energy Ep of each 40 K cloud. Error 
bars are statistical, with uncertainty smaller than the sizes 
of symbols for lower temperature data. Absorption images 
are shown for ksT/Ep — 0.35 (a) and 0.95 (b), including a 
black circle indicating the Fermi energy Ef. (c) A closer look 
at the fermion cloud shape reveals that it does not follow a 
Boltzmann distribution. The fit residuals of a radially aver- 
aged cloud profile show a strong systematic deviation when 
assuming Boltzmann (blue circles) instead of Fermi (red dia- 
monds) statistics. A degenerate Fermi cloud is flatter at its 
centre than a Boltzmann distribution, and falls more sharply 
to zero near its edge. 
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FIG. 4: Cross-species thermalisation. (a) The ratio of 
the temperature of 40 K to the temperature of 87 Rb approaches 
unity as 40 K temperature is lowered during sympathetic cool- 
ing, (b) We measure cross-thermalisation by abruptly reduc- 
ing the temperature of 87 Rb and watching the temperature 
of 40 K relax versus time. The data shown has an asymp- 
totic 40 K temperature of 27 /iK. (c) Measurements of OKFtb 
(diamonds) are compared to the "naive" model (dashed) and 
an effective range model (solid), both described in the text. 
The error bars shown in all parts of the figure are statistical, 
with the exception of temperature error bars in (c), which 
show the spread between initial and final 40 K temperature. 
The highest temperature point (open diamond) did not com- 
pletely thermalise, and is discussed further in the text. For 
reference, the s-wave crrtbRb is also shown (dotted). 



monicity at high temperature. Overall, both data and 
theory show that the 40 K- 87 Rb cross-section is reduced 
well below the 87 Rb- 87 Rb cross-section for a large range 
of temperatures, explaining the requirement for a slow 
initial RF frequency sweep for sympathetic cooling. Be- 
low 20 /iK, where no temperature lag is observed, <7KRb 
exceeds o- Rb Rb- 

We attribute the observed reduction in scattering 
cross-section to the onset of the Ramsauer-Townsend 
effect, in which the s-wave scattering phase and cross- 
section approach zero for a particular value of rela- 
tive energies between particles [27|. At higher tem- 
peratures, the scattering cross-section should increase 
again, however free evaporation from our trap limits 
our measurements to below 300 /iK. Additional partial 
waves may also affect scattering above the p-wave thresh- 
old of 110 /iK. Despite the high-temperature reduction 
in cross-section, 40 K and 87 Rb remain relatively good 



sympathetic cooling partners. Recent measurements of 
87 Rb- 6 Li sympathetic cooling Q, for instance, suggest a 
zero-temperature cross-section approximately 100 times 
smaller than <7KRb, i-e., a maximum cross-section roughly 
equal to the lowest value we measure here. 

The high collision rates in mixtures trapped with a 
/iEM enable us to cool fermions sympathetically to quan- 
tum degeneracy in 6 s, faster than previously possible. 
For comparison, direct evaporation of optically trapped 
6 Li to Fermi degeneracy has been achieved in 3.5s [28| . 
close to a Feshbach resonance. However, the all-optical 
approach does not lend itself easily to sympathetic cool- 
ing. The lower efficiency of direct evaporation requires 
a stronger source of 40 K, previously only achieved with 
multiple vacuum chambers 0, due to low available iso- 
tope abundance. 

In conclusion, we have achieved simultaneous quantum 
degeneracy of bosonic and fermionic atoms in a /iEM trap 
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and demonstrated an approach that can simplify future 
research with cold fermions. One prospect is the obser- 
vation of Pauli blocking in light scattering off degener- 
ate fermions 

mug. 

The high /iEM trap frequencies 
boost the ratio of Fermi energy Ep to the recoil energy 
H 2 k 2 /2rriK to ~2.5, within the range necessary to explore 
such quantum optical effects. 

Methods 

Loading. Our experimental cycle is similar to that 
described in Ref. with several key modifications 

emphasized here and in the main text. Approximately 
600 mW of incoherent 405 nm light desorbs 87 Rb and 
40 K atoms from the pyrex vacuum cell walls, boosting 
the MOT atom number one hundred-fold compared to 
loading from the background vapour. Potassium alone 
is first loaded into the MOT for 25 s, after which 87 Rb 
is loaded for an additional 3-5 s, while maintaining the 
40 K population. Both MOTs operate with a detuning 
of —26 MHz, until the last 10ms, when 40 K is com- 
pressed with a —5 MHz detuning. After MOT loading, 
3 ms of optical molasses cooling is applied to the 87 Rb 
atoms, while the 40 K atoms are optically pumped into 
the \F = 9/2, mp = 9/2) hyperfine ground state. 

Micro-electromagnet trap. 7-/xm-thick gold wires 
are patterned lithographically and electroplated on a sil- 
icon substrate. Two defects are present near the centre 
of the principal Z-wire, which result in the formation of 
three "dimples" in the trapping potential. We use the 
magnetic gradient generated by 30 mA of current through 
the U-wire to centre the magnetic trap on one of these 
dimples. 

Fitting absorption image data. Degenerate Fermi 
clouds are fit using a semi-classical expression for the 
optical density: Af2(Z exp [— g 2 /2r 2 ]), where g is the 
radial coordinate, Aji(Z) is the peak optical density, 
Z the fugacity, and f v (q) = — YliLi(~ l) is the 
Fermi-Dirac function. The temperature is given by 
k^T = r 2 mK/(wJ 2 + t 2 ), where r is the fit width and 
t the time of flight. Atom number is extracted using 
N = 2irr 2 f^(Z)A/ax, where o\ is the resonant absorp- 
tion cross-section. T/T-p can be extracted directly from 
the fugacity using (T/Tp)~ 3 = 6/3(2). Non-degenerate 
clouds are fit to a Gaussian distribution A exp [— g 2 /2r 2 ], 



with the same interpretation of r. Probes along both x 
and y (see Fig. ^| were used for imaging. Comparison of 
temperature measurements along axes of expansion sug- 
gest a 20 nK kick (possibly magnetic) is given to clouds 
along z, and that other temperatures agree systemati- 
cally at the 5% level. Data for residuals shown in Fig.|3J; 
is radially averaged about an ellipse defined by the two 
trap frequencies of the image plane. This ID radial data 
set is binned into 2-pixel bins, and fit as described. 

Scattering theories. The "naive" interaction model 
discussed in the text gives CKRb — 47ra 2 /(l + a 2 k 2 ), 
where a is the s-wave scattering length and k is the rel- 
ative wave vector in the centre of mass frame. Fig. 
shows the thermally averaged theory curves. Including 
the next order correction in the s-wave scattering am- 
plitude f(k) — — [1/a + ik + k 2 r c /2 + . . requires 
an effective range, which we calculate using 1261 to be 
r c = 20.2 ± 0.3 nm, for a K Rb = -10.8 ± 0.3 nm^O]. 

Analysis of thermalisation data. When the 87 Rb 
atom number Ajtb is much larger than the 40 K atom 
number, the relaxation of the 40 K temperature T to T^b 
is described by u = — uT~ 1 (l + mR\ J u/(rnR\ J + 'iTiK))2 (1 + 
w/2)~2 ; where u = (T/T^b) — 1, and thermalisation time 
r given by 

1 V2 CT K Rb y / m^m 2 ih uj 2 ± uji 

T on KB-fRb (mK + m-Rb) 

in which trap frequencies are for 87 Rb [2^. Fitting for r 
allows us to extract <XKRb- Note that all thermalisation 
data is taken with Nk below 4% of Ar^- 

Data in Fig. is analysed assuming a temperature- 
independent cross-section within the range of initial to 
final temperature. To check this assumption, we re- 
analyse the data using a self-consistent method that as- 
sumes an effective range temperature-dependence, and 
find a small upward shift of the best fit cross-section val- 
ues. Using this shift as an estimate of the methodology- 
dependent systematic error, we fit our four lowest tem- 
perature measurements with the effective range model, 
and find axRb — —9.9 ± 1.4 ± 2.2 nm, in agreement with 
Ref. [2(J. The second uncertainty reported is system- 
atic, and also includes uncertainty in the 87 Rb number 
calibration. 
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